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Abstract 



Explicit computations of tlie partition function and correlation functions of Wilson and 
Polyakov loop operators in theta-sectors of two dimensional Yang-Mills theory on the line 
cylinder and torus are presented. Several observations about the correspondence of two di- 
mensional Yang-Mills theory with unitary matrix quantum mechanics are presented. The 
incorporation of the theta-angle which characterizes the states of two dimensional adjoint 
QCD is discussed. 
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1 Introduction 



Two dimensional (2D) Yang-Mills theory is an interesting example of a topological field the- 
ory JI], H, ^. It has been used extensively as a toy model where features that it shares with 
higher dimensional gauge theories can be studied in a simplified context ^ • As a quantum 
field theory, it has no propagating degrees of freedom. It nevertheless exhibits a confining 
potential for heavy quarks at the tree level. It is also an example of a gauge theory where 
(in a planar geometry) the Migdal-Makeenko equations for Wilson loops can be solved ex- 
plicitly. Correlators of Wilson loops can be computed in closed form |^ and the area law 
associated with confinement physics can be demonstrated explicitly. Moreover, it has been 
demonstrated that the strong coupling, large limit of 2D Yang-Mills theory on the sphere 
can be rewritten as a random surface model p. This supports the long-standing conjecture 
that the infrared limit of a confining gauge theory is related to string theory. Its dynamics 
are also known to be intimately related with D=l unitary matrix models which are related 



to two dimensional string theories [|10| |ri|, ^ 



As a topological field theory, 2D Yang-Mills theory gives a field theoretical example of 
localization formulae 0]. It is a quantum field theory whose partition function is given 
exactly by evaluating the Euclidean functional integral in the saddle point approximation. 
The technique which is used to compute the functional integral is the so-called diagonalization 
method where gauge invariance is used to diagonalize the matrix- valued degrees of freedom to 
yield a simple, solvable model for the eigenvalues. The spectrum which is obtained this way 
differs from the conventional one. This difference is related to the existence of inequivalent 
quantizations of the gauge theory, depending on the order in which the constraints are imposed 
and quantization is done |]13|, |T^, [1^, 0], and was noted by Witten in his analysis of 2D Yang- 
Mills theory The diagonalization procedure has also been used to compute the partition 
function and correlators of Polyakov loop operators in ref. [0, 0. 

In this Paper we shall summarize some of our observations about 2D Yang-Mills theory. 
We shall begin by making a few remarks about the correspondence between Yang-Mills theory 
and unitary matrix quantum mechanics. We shall make extensive use of one such well-known 
correspondence 



19, 10 



20] , pT| where the unitary matrices are the Wilson loops. We shall also 
formulate another one where the dynamical variables are unitary matrices related to Polyakov 
loops. These matrices determine the gauge group holonomy that the wave-functions of heavy 
quarks accumulate in the quantized Yang-Mills theory. The latter correspondence is related 
to recent work where it was shown that the question of whether or not a gauge theory with 
adjoint matter and at finite temperature and density exists in a confined or de-confined phase 
could be related to the question of symmetry breaking in a certain nonlinear sigma model 
1T7| , In that model, the dynamical variable is a unitary matrix whose trace coincides 

with the Polyakov loop operator. The Polyakov loop operator is an order parameter for 
breaking the discrete symmetry and its expectation value characterizes confinement in finite 
temperature adjoint gauge theories [^, In general the action of the non-linear sigma 
model in question is a complicated effective action coming from integrating gauge and matter 
fields from the functional integral 



Tl that the 



However, it was argued in ref. 
effective action has a simple form in 2D Yang-Mills theory and the resulting model could be 
solved exactly using the diagonalization method [Q. This idea was generalized to Yang-Mills 
theory with external sources in [RTI] where it was used to show that there is a de-confining 



phase transition at infinite A^ limit of a hot and dense gas of heavy quarks. 
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In this Paper, we shall complete the program begun in ref. of presenting a full compu- 
tation of the correlators of loop operators in 2D Yang-Mills theory using the diagonalization 
method. We are particularly interested in the sectors of the theory with non-trivial theta 
vacua. We will present an explicit computation of the two-point correlators of Wilson and 
Polyakov loop operators in a theta-sector of 2D Yang-Mills theory at finite temperature. We 
shall consider both the case where the space is a circle, so that the spacetime is a torus, or 
the space is a line with fixed boundary conditions, so that the spacetime is a cylinder. We 
shall compute the partition function, find the behavior of loop correlators and deduce the 
quark- antiquark potential for both of these cases. We show that, in the limit of infinite vol- 
ume, our result reproduces the already known facts [^, that, in a sector of the theory 
with non-zero theta angle, fundamental representation quarks have a repulsive interaction 
and adjoint quarks have a screened, short-ranged interaction. Our results for the correlators 
of Polyakov loop operators are different from those found in previous computations of these 
quantities using the diagonalization technique |T^. 

Our principal result is the demonstration that the diagonalization technique can be mod- 
ified in such a way that it obtains the conventional spectrum of 2D Yang-Mills theory on a 
circle. Furthermore, we demonstrate that the two-point correlators which are obtained by 
this method agree in various limits with those obtained by the group theoretical character 
expansion technique. The distinct advantage of the diagonalization method is that it can be 
used to obtain explicit formulae for correlators of higher winding loops. 



1.1 Theta states 



It was first argued by Witten that 2D Yang-Mills theory, or 2D QCD with adjoint 



matter, has theta states. Since the Yang-Mills fields transform in the adjoint representation 
of the gauge group G, the true gauge group is the factor group of G by its center, C, G/C. 
Generally, the resulting group is not simply connected. For example, if G itself is simply 
connected {Jli{G) = 0). which we shall henceforth assume, then 

U^iG/C)=UoiC)=C (1.1) 

This has non-trivial consequences for any gauge theory that is defined on a spacetime which 
itself has non-trivial fundamental group. For the moment, let us assume that the space is a 
circle, 5*^.0 The group of time- independent gauge transformations, Q , which is formed from 
the smooth mappings of the circle. Si, to G/C has the property 

Uo{g) = MG/c) = c (1.2) 

Since the Hamiltonian is gauge invariant, and the quantum states are invariant under gauge 
transformations in the component of the gauge group Q which contains the identity, the 
quantum states must carry an irreducible unitary representation of the Abelian group C. All 
such representations are one-dimensional. Each distinct representation is called a theta-state 
(this terminology derives from the representations of elements of U(l) by phases e*^). 

^ We could as well use the real line R^. In the case of the real line, the gauge group is taken as the set of 
those gauge transformations which go to the identity at infinity. 
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All gauge invariant operators have vanishing matrix elements between states which carry 
different representations of C. Thus, there is a superselection rule which allows one to choose 
one particular representation to characterize all of the quantum states of the theory. 

For example, if G = SU{N) then its center is C = Z^, the cyclic group of order N and 

Uo{g) = Hi {SU{N)/Zn) = (1.3) 
There are N irreducible unitary representations of the center, Zat, which have elements 

1, e^'\ . . . , e'^^~^^^ 
where the N allowed values of 6 are 

e = 0,2TT/N,...,2nk/N,...,2n{N -1)/N 



The partition function of Yang- Mills theory is 

Z = jdA^ exp (- j^^ d\^TT{FlS) (1-4) 

Let us consider the case where the Euclidean space-time is a Riemann surface, M. The 
gauge fields belong to topological sectors. In a generic topological sector, the connection one- 
form A is not globally defined, but is a smooth 1-form only in coordinate patches. Consider 
a triangulation of the space M by a set of simply connected subspaces labeled by I 
(M = UiT^) and let the connection in patch / be . At the boundary between patches I 
and J, the connections are related by a gauge transformation A^ = {A-^)^''' . 

We can change a gauge field which is in one topological sector to one in another sector by 
the following surgical procedure Consider a particular field configuration Aj in patch / 
and transition functions gjj which relates Aj to the connections in neighboring patches. We 
can change topological sectors by replacing Aj by its gauge transform A^ and the transition 
functions gu by gjjg if we choose g to have the property that the restriction of g to the 
boundary of the patch (which is isomorphic to S^) is an element of Q, the smooth mappings of 
S*^ to G/C, which is in a non-trivial component of no(^) = C. The resulting field configuration 
will be in a different topological sector from the original one. The number of topological 
sectors is equal to the number of elements in C. Summing over the topological sectors in 
the path integral with complex unimodular weights given by the corresponding elements of 
the irreducible unitary representation of the center, C, produces the partition function in a 6' 
sector. 

The above arguments depend on certain smoothness requirements for the class of field 
configurations which contribute to the path integral and also of the gauge transformations 
which are allowed in the Hamiltonian formalism. However, as independent evidence for the 
existence of theta states, they have been given a physical interpretation similar to the theta- 
vacuum in the Schwinger-model []35| , |36| as being created by the existence of charges at the 



boundaries of the space. The theta parameter is known to affect the spectrum of mesons in 
2D QCD with heavy quarks ||3^, |3^. Arguments of stability of such systems have also been 
shown to be consistent with the topological classification of theta states P^, 53. 
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1.2 Overview and Summary of Results: 

In Section 2, we shall review the quantization of Yang-Mills theory in 1+1- dimensions when 
the space is a circle. We pay particular attention to the existence of large gauge transfor- 
mations and the associated appearance of 9 vacua. Then, in Section 3, we derive a number 
of different representations of the heat kernel of 2D Yang-Mills theory. We also present an 
explicit calculation of the heat kernel using the diagonalization technique. Section 4 contains 
a computation of loop correlators on the cylinder and torus. We also obtain results for SU(2) 
and compare them with a similar computation using character expansion techniques and find 
that they agree. 

We summarize our results as follows: 



We obtain the partition functions and the correlators of two loop operators for 2D QCD 
on both the cylinder and the torus in a 6'-state when G = SU{N) given in formulae 
( p.36| ), ( p.53| ), ( |4.8| ) and ( [4.12| ). This is done using the representation of the heat kernel 
in (|1.5|) below. In the process, we show that the trace of (|1.5| ) leads to the character 
representation of the partition function ( |1.6| ). 



It is well known |[T9[1 that the physical states of 2D Yang- Mills theory on a circle {x G [0,L]) 
are class functions of the Wilson loop 



where f f ^ = 1 and 



^phys(^) = tpphysivUv^) 



U{t) = Pexp J dxAi{x,T)^ 



At finite temperature where the Euclidean time is periodic {t G [0,r = 1/T]), the Polyakov 
loop operator is also of interest. 



g{x) = Pexp I z J dtAo{x,t) j 



The heat kernel of 2D Yang-Mills theory, defined by 



{^2 



-Ht 



[dU^][dU2]r2iU2)MUl)K[T■,U2,U^ 



so that 



K[t;U2,Ui] = (f/2 



-Ht 



is defined as a symmetric function of the eigenvalues of U2 and Ui separately. It can be written 
in three ways: It has the well known decomposition in terms of characters Xr{U) = Tr^jf/ of 
U in irreducible representations R, 

K[t;U2,U,] = Ex^j(f/2)Xii(f/i)e-^'^^^^(^) 

R 

where C2{R) is the second Casimir invariant of R. In addition, we shall find that it can be 
presented 
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• as a D=l unitary matrix model with a background gauge field, 

K[t; U2, Ui] = J [dg{x)] exp (^"^ £ dxTi \Vg + iA2g - igAi\^'^ 
where (?(0) = g{L) and 

D=l gauge invariant matrix model 

K[t;U2,Ui] = J[dU{t)]dA{t)exp(^-^J^ dtTT\u -i[A,U]\^y 

■ 5ci(f/(r), f/2)5ci(t/(0), U,)J{U{t))J{U{0)) (1.5) 

where J{U) is the Vandermonde determinant and 6ci{U, V) is the conjugation invariant 
delta function which equates the spectra of U and V. 

• The partition function of Yang-Mills theory in a 6'-state on the 2-sphere with area Lt 
can be written as the D=l unitary matrix model 

Z[sphere;^] = J2Z{z,e) J[dg,]exp (^-^ J^ilVgf^ 

where ^'2(0) = zg^iL), z is an element of C, the center of G, and Z{z, 9) is the element 
of the 6'-representation of C corresponding to z. 

• The partition function of Yang-Mills theory in a 6'-state on the 2-torus with area Lt 
can be written as the gauge invariant D=l matrix model 

Z[torus;^^] = ^Z(^,^)|[rf(72][rfA]exp(^-^|^\r|V(7.-2[A^7j|') 

In the equivalent character representation the partition function of Yang-Mills theory 
on a Riemann surface of genus g and in a ^-sector is given by 

Z[g^ e]=Y, (dim R)^"^^ 5{R, Re) exp (-eVLC2(i?)) (1.6) 

R 

where Rq are those representations which have the property 

xrA^u) = z{z,e)xRAU) (1.7) 

2 Hamiltonian formulation of Yang-Mills theory on the 
circle 

In this Section we shall review the Hamiltonian formulation of 2D Yang-Mills theory. Canon- 
ical quantization of the action^ 

S = j dxdt ^Tr (ii - VAof (2.1) 



Note that we use an unconventional normalization of the charge 
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yields dynamical variables which are the spatial component of the gauge field Ai{x) = A{x) 
and the electric field E{x) = A{x) — Vv4o(x), both of which are matrices in the fundamental 
representation of the Lie algebra of the gauge group, and can be expanded in a canonical set 
of generators as 

A{x) = ^TM"(a;) , E{x) =Y,T''E°-{x) 



For a compact semi-simple Lie algebra, 



1 



~ 2 



^ j!abcrpc 



The Hamiltonian is ^ 

H = 2e^J^ dxTi(^E^{x) 

and the non-vanishing canonical commutation relation is 

'A''{x),E\y)] =i6''^6{x-y) 

Aq is a Lagrange multiplier which enforces Gauss' law as the constraint is 

VE{x) + i[A{x),E{x)] ~ 

and all variables have periodic boundary conditions, 

E{L) = E{0) , A{L) = A{0) 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



The gauge transformation is 

E{x) E3{x) = g{x)E{x)g\x) , A{x) A^{x) = g{x) {A{x) - iV) g\x) (2.6) 

The periodic boundary conditions for A{x) and E{x) are preserved by gauge transformations 
where the matrix g{x) and its first derivative are periodic up to an element of the center of 
the gauge group, 

g,{L) = zg,{0) (2.7) 

where z & C. When z = 1, gi{x) is said to implement a "small gauge transformation" and 
when z ^ 1, gz{x) is said to implement a "large gauge transformation" . 

The operator on the right-hand-side of the constraint ( p.4|) generates infinitesimal small 
gauge transformations. As a result of the constraint ( p.4| ), the physical phase space is the set of 
equivalence classes of configurations E{x) and A{x) where fields are in the same equivalence 
class if they are related by a small gauge transformation (i.e. one with a periodic gauge 
function and with z = 1 in ( p.7|) ). The constraint could be solved at the classical level 
by choosing representatives of the equivalence classes. This procedure has been discussed in 
|]37| , |38| . Alternatively, after finding the quantum realization of the un-constrained theory with 
commutator (|2.3|) and Hamiltonian (|2.2| ), the constraint (|2.4|) could be imposed as a physical 
state condition which chooses a subspace of the quantum states as 'physical states'. It implies 
that physical states of the theory are invariant under all gauge transformations which can be 
generated by iterating infinitesimal transformations, i.e. all small gauge transformations. The 
coset group of all gauge transformations modulo periodic ones is isomorphic to the center 
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of the gauge group, C. The physical states must transform under an irreducible unitary 
representation of this coset. 

In the functional Schrodinger picture, where states are wave-functionals, ^E'[A], of gauge 
field configurations and the electric field is realized as 

the physical state condition, 

{VE{x) + i[A{x),E{x)]) ^phys.[v4;0] = , (2.9) 
implies that the wave-functionals of physical states transform as 

^^^,ys.[gz{A-iV)gl-,e] = Z{z,9)^^^,,\A-e] (2.10) 

where gz{x) has the boundary condition in ( p.7| ) and Z{z^ 9) is the representative of the center 
element z in the representation labeled by 9. Since the center of the gauge group is an Abelian 
discrete group, all irreducible unitary representations are one- dimensional and the number of 
inequivalent representations is equal to the order of the group. The solution of ( |2.1(]| ) is 



^phys.[A;^^] = i^[U;9] (2.11) 

where U is the unitary matrix 

U = V exp (^i £ dxA{x)j (2.12) 

with the additional constraint that the wave-function is a class function of U : 

i/j[U;9] = ^[hUh^;9] (2.13) 

for any h in the fundamental representation of G. This implies that the wave-function is 
a symmetric function of the eigenvalues of U. On the wavef unctions, we must impose the 
further requirement that they lie in a particular theta-sector, 

iP[zU;9] = Z{z,9) iIj[U,9] (2.14) 

Let us consider the operator Un constructed from the path-ordered phase integral with the 
gauge field in an arbitrary irreducible representation R of the Lie algebra, 

Ur = Vexp(^i £ dxA''{x)T^ (2.15) 

where are the generators in representation R. The electric field operates on Ur as 

E'^ix) U = Vexp £ dwA\w)T^ V exp {i £ dwA%^[w)^ (2.16) 

and the field space Laplacian which is in the kinetic term in the Hamiltonian as 

Ur = C2{R)Ur (2.17) 
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SA^ix) 



where C2{R) is the second Casimir invariant corresponding to the representation R which is 
obtained from the formula 

C2{R)-X = j:nTR (2.18) 

a 

We thus find physical eigenstates of the Hamiltonian by forming traces of Ur, i.e., the char- 
acters which are defined by 

XRiU) = Tr Ur (2.19) 

so that 

H XRiU) = e^L C2{R) Xr{U) (2.20) 

Thus, the action of the Hamiltonian on the characters is proportional to that of the group 
Laplacian, A(G) whose spectrum is the set of second Casimir invariants of irreducible repre- 
sentations. 

Characters are orthonormal with respect to integration over the Haar measure, 

J[dU]xR{U)xR'{U) = 5r,r, (2.21) 

This inner product can be obtained from the natural inner product, which is functional 
integration over A, by gauge fixing, 

J ciA(x)*;hys.[^]*phys.[A] = const. • J [dU] r{U, e)i;{U, 6) (2.22) 

Products of characters are also important as they effectively carry out the multiplication 
of irreducible representations 

XrAU)xrAU) ^ Xr^^rAU) (2.23) 

If i?i ® i?2 = ®N^tR,R3 

is a decomposition into irreducible representations of the product 
then by linearity of characters 

XR, {U)XR. {U) = m^R^XRs (U) (2.24) 
Consequently we have the definition of the fusion numbers, 

= / mxRAU)XR.{U)XR,{U) (2.25) 

Thus the set of all states is identified with the characters of U. The theta-states in a 
given theta-sector are a subset of these. For a certain representation of the center of the 
gauge group, with elements Z{z, 9) we must chose the characters for those representations of 
G which have the property 

Xr{zU) = Z{z,e)xR{U) (2.26) 

The spectrum of the Hamiltonian is then the set of second Casimir invariants corresponding 
to the representations R of the group G whose characters have this property. 

In the following we shall deal mostly with the special case where the gauge group is SU(N). 
Then the A^^ — 1 generators T" are NxN traceless Hermitean matrices. The center of SU(N) is 
the cyclic group of order A^, Zn. The representations which have the property Xij(e^'^*"/^t/) — 



8 



Z{6)xr{U) = e*"^X-R(^) with 9 = 2TTk/N are those representations whose Young tableaux 
have k modulo N boxes. If we consider U{N) representations, this is equivalent to keeping 
only those with linear Casimir Ci{R) = kmodN. 

It is well known that 2D Yang-Mills theory is equivalent to a particular version of unitary 



matrix quantum mechanics. This can either be seen by fixing a particular gauge as in [10, 11 



or by the following argument which shows that the Hilbert space of states, energy levels 
and degeneracies of the two theories are identical of we identify the unitary matrices as the 
Wilson loop operator in ( [2.121 ). The Hamiltonian of unitary matrix quantum mechanics is 
proportional to the group Laplacian, 

Hqm = e^LA{G) (2.27) 
can be derived by canonical quantization of the action 

S = ^JdtTr (f/tf/) = ^JdtTr (zT"t/t(t)f7(t))^ (2.28) 

The canonical momenta are given by 

n° = ^4^Tr {iT^U^u) (2.29) 

and the Hamiltonian is 

Hqm = e'Lj2m' = e'LAiG) (2.30) 

a 

The canonical momentum operators have the Lie algebra 



ir'^'w (2.31) 



and 

[n^[/] = T^U , [n",f/1'] = -U^T" , (2.32) 

The wave-functions of this system are the group elements in unitary irreducible represen- 
tations of the group, 

A{G) Un = C^iR) Ur (2.33) 

The inner product is given by the integration over the group with the invariant Haar measure, 
where the wave-functions have the property 

J [dU] {Ul)ki {Ur\, = -^-^ Sa5,k (2.34) 

Each component of the unitary matrix in the representation i? is a linearly independent, 
normalizable wave-function. The degeneracy of each eigenstate is equal to the number of 
linearly independent components, i.e. (dim i?)^. 

However, to produce the spectrum and degeneracies of 2D Yang-Mills theory, the wave- 
functions must be restricted to class functions of group elements so that, for each representa- 
tion of G. only the wave-function Xr{U) is allowed. At the operator level, this restriction is 
realized by the constraint 

Y^u" ■ tt (t'^t'' - rm'^u^) ~ o , (2.35) 
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The operator in this constraint generates infinitesimal adjoint transformations of the group 
elements. This constraint can be enforced by a Lagrange multiplier and the Hamiltonian and 
constraint can be obtained by canonical quantization of the action 

S = ^ldm\u{t)-i[U{t),Ait)]\' (2.36) 

This action has a gauge invariance under 

A-.A^ = gAg^-tgg^ , U-^W = gUg^ (2.37) 

Here, the gauge field A{t) is a Lagrange multiplier enforcing the matrix mechanics analog 
of Gauss' law. The constraint implies that the adjoint action of the group generators on 
the wave-function vanishes. Wave-functions which are annihilated by the constraint are class 
functions of the group variables, 

m] = ^[gUg^] (2.38) 

Of the (dim(i?))^ linearly independent eigenfunctions Ur of A(G), only one, the character 
Xr{U) = TtUr has this property. 

In this theory, the expectation value of the Wilson loop operator in a given quantum state 
is given by the fusion number, 

(i?i|x«(f/)|i?2) = l[dU]x*n,{U)xR{U)xR,iU) = N^^j, (2.39) 

Similarly, the correlator of a product of Wilson loop operators in arbitrary representations is 
given by 

{R\XRAU)XR,{U)...XRdU)\R') = J[dU]xR{U)XR,iU)XR,{U)...XRdU)XR'{U) (2.40) 

For the group SU(N), explicit formulae for these moments can be obtained. Also, the partition 
function for this theory at finite temperature T = is given the the familiar expression for 
the partition function of Yang-Mills theory on the torus, 

Z[t] = exp (-e2LrC2(i?)) (2.41) 

R 

We would expect that this partition function could also be written as a functional integral 
for the d=l unitary matrix model with action (p.36 ). We shall show in the following sections 
that, modulo some subtleties with boundary conditions, this is indeed the case. 



3 The heat kernel: Yang-Mills theory on the cylinder 

In this Section, we shall give two equivalent presentations of the heat kernel for 2D Yang- Mills 
theory. The first shows the connection between 2D Yang-Mills theory and a certain principal 
chiral model which generalizes a result in [T^. In this representation, the dynamical variables 
are the unitary matrices which would parallel transport heavy quark wave functions from the 
initial to the final (Euclidean) time. If we consider a partition function, so that Euclidean 
time is periodic, their traces would be Polyakov loop operators. 
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Alternatively, we shall find it useful to represent the heat kernel as a c = 1 unitary matrix 
model where the matrices are Wilson loop variables. We shall show, with details in subsection 



0| , the equivalence of the latter model with the well known character expansion of the heat 
kernel. 



We begin by considering the propagation function in Yang-Mills theory. 



K[t;A2,A,] 



Ai 



(3.1) 



where \A) is an eigenstate of the gauge field operator A"-{x), H is the Hamiltonian in ( |2.2D 
and P is a projection operator onto gauge invariant states. For the moment, we consider 
states which are invariant under periodic gauge transformations only. The projection can be 
implemented by gauge transforming the field Ai at one side of the trace and integrating over 
all gauge transformations. 



K[t;A2,A,] = l[dg{x)]{A2\e-''-\Al 



(3.2) 



(where we normalize the measure so that j'[dg{x)] = 1 and g{0) = g{L).) The integrand in 
( |3.2| ) is the heat kernel which obeys the following equation. 



(a;). 



A, 



-Ht 



Ai 



with the boundary condition 



lim ( A2 



-Ht 



Al) = l[5iA2{ 



x] 



Af{x)) 



(3.3) 



(3.4) 



These are solved byQ 

K[t; A2, A,] = J [dg{x)] exp (^-^ dx Ti^A^ - A^f^ 
We can re-arrange the action in ( p.5|) to put it in the following form, 

K[t- A2, Ai] = J [dg{x)] exp dxTi \Vg + ig{x)A^{x) - iA2{x)g{x) 



(3.5) 



(3.6) 



This is the path integral for a 0+1-dimensional principal chiral model with external gauge 
fields where we treat the spatial variable x as Euclidean time. With this identification, 
it also coincides with the partition function unitary matrix quantum mechanics coupled to 
external gauge fields. Since the Haar measure has the properties [c?(7(a;)] = [d{u{x)g{x))] = 

^Here, we have dropped a zero point energy term for the Hamiltonian. Also, this equation should be 
divided by the normalization factor {2T:e'^T)^^^^\ If we use zeta-function regularization, this normalization 
has a very simple form. Since 



LS{0) = linr (1 + 2 ^ 1/n') = 



the normalization factor is one. 
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[d{g{x)v'^{x))] where u{x) and v{x) are unitary matrices, the heat kernel has the property that 
it is invariant, under gauge transformation of Ai{x) and ^2(0;) separately, 

K[t;AIA^,] = K[t;A2,A,] (3.7) 

This implies that the heat kernel is a class function of each of the two unitary matrices 

Ui = V exp (^i dxAi{x)j , U2 = V exp (^i £ dxA2{x)^ (3.8) 

Consequently we define 

K[t; A2, A,] = K [t- U2, f/i] = {U2 |e-^^P| U,) (3.9) 
which has the invariance property 

K[t; U2, f/i] = K[t- uU2u\ vUiV^] , (3.10) 



As is shown in the subsection p.2| , this form of the heat kernel (|3.9| ) has a natural expansion 
in terms of group characters, 

ir[r;f/2,f/i] = Xr{U,) exp(-eVLC2(i?)) Xr^ (3.11) 

R 

which, with the expression ( |3.8| ) and ( |2.33| ) can be seen to satisfy the heat equation 

(J- + e^LA{G)^ K[t; U2, U,] = (3.12) 

where the gauge group Laplacian A(G) operates on t/2. Also, the boundary condition, 

K[0,U2,Ui] = [[dg{x)]Y[6{A2{x)~Al{x)) = (const.) ■ (5d(t/i, t/2) (3.13) 

The delta function on the right-hand-side of (p.l3|) is the conjugation invariant delta func- 
tion which equates the eigenvalues of the two unitary matrices. It can be defined by group 
integration, 

Sci{9i,92) = f [dV]6{VgiV^g2,I) = Y.Xr{9i)x*r{92) (3-14) 

R 

Furthermore, the heat kernel for 2-dimensional Yang- Mills theory on a cylinder of length r 
and base circle L can be written in terms of the partition function of a gauged principal chiral 
model with open boundaries. This, as we showed in Section 2, is equivalent to 2D QCD if 
we restrict our attention to Wilson loops. The path integral representation of the heat kernel 
for that theory should have a standard form using action (|2.36| ) and integration variables A 
and U. Some care must be taken to ensure that the heat kernel is a correct function of the 
eigenvalues of the Wilson loop operators. In particular the sewing property, 

Z[t; U2, f/i] = J dUiu)Z[u; U2, Uiu)]Z[T - u; U{u), U^] , (3.15) 

must be satisfied. 
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The conditions that fix the loops at the two boundaries of the cyhnder (say U (0) and U{t)) 
to be Ui and, up to an element of Z^, U2, can be imposed introducing the corresponding 
delta-functions: 

Z[r;Ui,U2] = f n [dA{t)][dU{t)]e-^^o''''^\^-'^^'''^\' 
te[o,T] 

5ei(t/(0),t/i)(5ci(f/(r),f/2)V^(f/(0))^(f/(r)) . (3.16) 

The factors ipiU^O)) and ipiU^r)) are boundary wave functions that have to be introduced in 
order to guarantee that the sewing prescription ( p. 151 ) for the kernel is satisfied. By means of 
eq.( |3.16|) , the left hand side of eq. (|3.15|) reads 

^ dU{u) n [dA{t)][dU{t)] n [rfA'(t)][rff/'(t)]e-i^/o"'^*Tr|c/-^[A,i/]|' 

te[0,M] ig[M,T] 

e-^L />*Tr|l/'-.[A',C/']|^^^^ ^^^Q^^ (^^^^^ ^^^^^ ^jj,^^^^ jj^^y^ 

Sci {u\t), U2) ij{ummiumu\u))ij{u'ir)) . (3.17) 

We will show in the next section that the integration over one of the A variables at the 
point u, produces the squared inverse of the Vandermonde determinant for a unitary matrix, 

J{U)= n2sini(r-/) (3.18) 

which only depends on the eigenvalues of f/, exp (20"), a = 1, . . . , A^. Consequently, Eq. (|3.15|) 
will be satisfied if we choose iIj{U) = J{U), so that 

K[t- U2, U^] = J [dA{t)][dU{t)] exp (^-^ £ dt Tr \u{t) - i [A{t), f/(t)]|') ■ 

■ 5ei(f/(0),f/i)4i(f/(r),f/2)J(f/(0))J(f/(r)) . (3.19) 

The formulae (|3.6|) and (|3.19| ) are two alternative expressions for the heat kernel, one in 
terms of matrix quantum mechanics with a background gauge field and the other in terms of 
gauge invariant matrix mechanics with the external matrices appearing as boundary condi- 
tions. 

The propagator in a 6'-sector can be obtained by projecting either the initial or final state 
of the heat kernel onto ^-states. This is done by summing over all transformations of f/i or 
U2 by elements of the center of G, C. and, in that sum, weighting each term by the phases 
Z{z, 9) corresponding to the 6'-representation of C, 

K[T,e-U2,Ui] = ^Z(^,^)ir[r;zf/2,[/i] (3.20) 

For the three presentations of the heat kernel given above, this has the effect 

K[t, 6; U2, U,] = Y.6 {R, Re) X^^mXR^ ■ exp {-e'LTC2{R)) (3.21) 

R 

= Y^Z{z,e) J[dg^{x)] ^^p(^~^^_^ dx TT\Vg^{x) +ig^{x)Ai{x) -iA2{x)g^{x)\^^ 

(3.22) 

= J2Z{z,e) J[dAit)][dU{t)] expi^~^£dtTT\u{t)-z[A{t),U{t)] 

■ 5ci(f/(0), Ui)6,,{U{t), zU2)J{U{0))J{U{t)) (3.23) 
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Here, Rg are those irreducible representations of G which have the property 



In the next section, by means of a diagonahzation technique, we shall show that the 
functional integral in (|3.23| ) can be actually performed so that the heat kernel can be expressed 



directly in terms of the eigenvalues of t/i and 1/2- 



3.1 Calculation of the heat kernel by diagonahzation technique 



In Ref. ||12|| two dimensional Yang-Mills theories were written in terms of a Kazakov-Migdal 
model and the heat kernel on the cylinder was computed using a diagonahzation procedure. 
We shall follow here an analogous procedure for the calculation of the partition function ( |3.23D 
of the gauged principal chiral model. In fact, the integral in (|3.23|) can be done by using 
gauge symmetry to diagonalize at each point the unitary matrix U(t) = V(t)U^(t)V^^{t) 
where U^{t) = diag(e*'^ W^g^p e**^ *^*^). The Lp variables are angles since G [0,2??]. 

The measure in the integral has the form 

mt)]= n i[dv''{t)j'{vm(Y.^''it)]mt)] , (3.24) 

tG[0,r] " \ a / 

where J{ip{t)) is the Vandermonde determinant for a unitary matrix, and the action is 

^ = ^ / ^ + ^ Mle^""" - ^^"'1' • (3-25) 

The integral over Aap, where a jS, cancels the Vandermonde determinant in the integration 
measure. The integration over the diagonal components of A yields an infinite factor which 
compensates the infinite normalization of the plane- wave states which were used in ( p.9|) . 
The kernel ( ^^231) becomes an integral over the eigenvalue variables. In the ^-sector using for 



J2z^iz^(^) the representation X^n 6xp(— m6'fe) with 9k = 2TTk/N, we get 

n=0 '' tg[0,T] a=l \a=l J 

J dV{0)S{V{0)e'^^°^V\0)e-'^\l) J dV{T)6 {V{T)e'^^^^V\T)e~'^^+^™/^ , l) 

(3.26) 

where e*^^'^ denote the diagonal forms of Ui^2 and the integration measure integrates if{t), 
at each point t, over the range [0,27r]. The integral is invariant under the field translation 
symmetry 

<^°(t) if'^it) + 27m° (3.27) 

where is an integer.0 This symmetry can be used to extend the limits on the integration 
over (p{t) at each t to infinity. 

^Actually in this case Ua can be any constant, but when computing loop correlators only the symmetry 
( |3.27 ) with integer Ua will survive. Note that this is not a symmetry of the action or the integration measure 



in the path integral (3.23) separately but only appears after diagonahzation. 
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The integrals on ^"'{t) with t G (0, r) can be performed by a ^-function regularization 
of the divergences. First one can integrate on (p^{t) taking advantage of the 5(Z]a V^"(^))- 
Changing the functional integration variables to the real and periodic (p°'{t), <^"(t) = <^°(0), 
by means of 



the integrals on the open interval t G (0, r), give 

N-l 



exp 



E (y."(r)-(p"(0))fi„,(/(r)-(/,'^(0)) 



where 



N-l 

n n^rwexp 

iG(0,T) 0=1 



" VOL G 
and n is the (A^ — 1) x (A^ — 1) matrix 



2e2Lr Jo 



N-l 

dt J2 ^ ^a/3<^ 
0/3=1 



(3.28) 



(3.29) 



(3.30) 



/2 


1 . 




1 


2 . 


. 1 


u 


1 . 


• 2) 



(3.31) 



Expanding V5"(t) in modes 



2 00 

fc=— 00 



a i2iTkt/T 



(3.32) 



the functional measure is defined as Y{tdip{t) = H^dak. The integration over the zero modes 
produce an infinite, irrelevant, temperature independent factor proportional to the volume 
of the gauge group. The functional integral in Za is proportional to the determinant of the 
Laplacian, 



Za=l[[ det 



271 



det 



e^Lr 
27r 



-C(o) 



,(iv-i)c'{o) 



. (3.33) 



where we have used zeta-function regularization with 

00 1 



k=l 



^2s 



(3.34) 



and C(0) = -1/2 and C'(0) = -(log27r)/2. 

We are now left with the integration on the (/9"(t), a = 1, . . . , at the boundaries t = 0, r. 
This can be easily done by rewriting the conjugate invariant delta functions in (|3.26| ) according 



to 12, 



dV5 (Ve'^V^e 



iX 



+00 (_l\P+iN-l)J2'^^^n^ N 



(3.35) 



0=1 
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where P denotes a permutation of the indices and (—1)^ its parity. By substituting Eqs.( p.2^^ ), 
(|3.33| ), ( p.35| ) into Eg. (|3.26| ) and integrating in V5"(0), v^"(r) we finally get 



K[T,ek;U2,U,]=NWN 



e^Lr 



^ N-l 



N 



5 ( ^ n^, + n I exp 

{na}=-a 



\a=l 



EE 

n=0 P 

N 



-2mn[k-N{N-l)/2]/N 



2e^LT 



a=l 



E (^A? - A?") + 27m„ + ^ j 



27rn\' 



(3.36) 



The delta-function imposing the constraint J2a na + n = arises from the S{J2a V^«) by taking 
into account that the U2 eigenvalues are twisted by an element of the center and that, being 
the gauge group SU (N) , we can choose J2a Aq = for the loops at the boundaries. We shall 
show in the nex subsection that Eq. (|3.36|) is equivalent to the character representation (p. 21 



3.2 From the functional integral to the character expansion 

Following Ref. ||12||, it is possible to show that, by means of a Poisson resummation, Eq. (|3.36|) 
can be rewritten in terms of a character expansion. It is convenient first to rewrite the delta 
function on the integers na as an integral 



M +^ = / — exp Ui* > ra„ + 



n 



(3.37) 



One can then complete the square in n^, so as to rewrite Eq. ( |3.36| ), according to 



K[T,9k;U,,U2] = N\ 



N-l 



EE 

71=0 P 



^_^^P^-2TTin[k-N{N-l)/2]/N 



J(Ai)J(A2) 



+00 

E 

{na}=-oo 



27T 



exp 



2.- 



e'^Lr 



27r 



(3.38) 



where we used that Yl,a=i A^ = 0, i = 1, 2. Applying the generalized Poisson resummation 
formula 



exp (—rrgijm^m^ — 27ri 



mm a, 



1 



m=— 00 



det gij 



exp —Tig'^ {rii — ai){nj — aj) . (3.39) 



where g'^^ is gij inverse, and then completing the square in Eq. (|3.38|) becomes 



K[r,d,-U^,U2] = (e'LrNf" E E E 

{m„} n=0 P,P' 

TV 



Af-1 



(_l)(-P+-P')e-2™[A;-Af(Af-l)/2]/Af 



J(Ai)J(A2) 



Em„ (a^") - Af") + 2vrml) 



exp 



exp 



e^Lr 



TV 



a=l 



N 



/ N 

E 



rrir 



Ne^Lr 



N 



87r2 



(3.40) 
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where we introduced a redundant sum over permutation cancehng the A^! factor. Note that 
Eq.( |3.40D is symmetric with respect to any permutation of the rria and vanishes when any 
two of rria, m/3 are equal. Therefore one can introduce an ordering between the integers 
mi > m2 > . . . > rriN multiplying by a factor A^!. One can further notice that all the terms, 
except the last integral are invariant if all the rria are shifted by the same integer. One can 
then define the integers Tq, = rria — — N, a = 1, . . . ,N and use the sum over to 
complete the Gaussian integral in by shifting ^ — 27Tmj\f — 2'kN . These manipulations 
lead to 



N-l 

K[T,ek;U,,U2]=N\ 

n=0 T-i>r2...>r]v = 



A^(A^- 1) 



1 dete*'""^^ dete" 



27ry 



J(Ai) J(A2) 



exp 




N 



(3.41) 



The ingredients of Eq.( 3.41| ) are now easy to recognize. The last exponent is, up to a constant, 
the SU{N) quadratic Casimir in the representation R whose Young table has Tq, + a boxes in 
row. 



the a*^ 



N 



N 



a=l \a=l 



2C2{R) + 



12 



(3.42) 



The second term in the r.h.s. of Eq.( |3.42| ) is the trace of the inverse Cartan matrix and is 
related to the scalar curvature of the group manifold. It provides the energy of the lowest 
representation. The first exponential in Eq.( |3.41| ) projects the sum to those representations 
Rk whose Young tableaux have rows with a number of boxes equal to the theta vacuum k, 
mod[iV]. Moreover the determinants are related to the character of the matrix ?7i 2 in the 
representation R. 



Xr{U) 



det (. 



J(A) 



(3.43) 



Up to an irrelevant, temperature independent constant, one can then write the final equation 



as 



K[t, Ou; t/i, f/2] = E XrAUMrM) exp 

Rk 



e^LT 



2C2{Rk) + 



12 



(3.44) 



3.3 2D Yang-Mills theory on the sphere as matrix quantum me- 
chanics 

Once one has the heat kernel for 2D Yang-Mills it is an easy matter to construct the partition 
function for the theory on a sphere or a torus. If one imagines shrinking the ends of the cylinder 
to points then the resulting topology is that of a sphere. In terms of the the discussion of 
the previous section this corresponds to setting Ui = U2 = I. Consequently, from ( p.ll|) the 
partition function on the sphere is 

ZIS"^; r, 6] = K[t; 6,1,1] (R^ Re) (dim(R))2 exp (-eVLC2(R)) (3.45) 

R 

where we have used the fact that x_r(1) = dim(R) and Rq is defined in (HI). 
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More interestingly, if we carry out the equivalent procedure with the path integral principal 
chiral model of the heat kernel (|3.6| ) and set Ai{x) = ^2(0;) = we find the following form 

Z[S';T,e] = J2Z{z,e) J[dg,{x)] exp(^-^ j^'' rfxTr |V(?,(x)|'^ (3.46) 

where g{L) = zg{0). 

Consequently we can interpret the partition function for the sphere as unitary matrix 
quantum mechanics. 



3.4 2D Yang-Mills theory on the torus as gauge invariant matrix 
quantum mechanics 

Similar to the case of the sphere it is easy to imagine manipulating a cylinder to form a torus 
- identifying the ends and sewing them together - which produces the trace of the heat kernel. 
This process is easy to carry out in the character representation by identifying Ui = U2 and 
integrating. Hence we have the partition function on a torus of area tL 

Z[T'^- T,6]= j [dU] K[t, e;U,U]= ^e) ^^P {-e'^TLC2{R)) (3.47) 

R 

Explicitly, for the special case of SU(N), where representations are denoted by the usual 
Young Tableau row variables h > h > ■ ■ ■ > > with I = we can give an explicit 
form for the partition function on the torus. When 9 = 2T[k/N, we consider only those 
representations where J2ih = k (modA^). Then 



Z[r,ek = 27Tk/N] = ^5Ar(^/„A;)exp 



N 

-e^TL + N + l-2j - l/N) 



(3.48) 



Equivalently, in the following section, we shall show how this result can be obtained by 
taking the trace of the heat kernel represented as a functional integral for the matrix model 

dm). 



3.5 Partition function on the torus in the functional integral ap- 
proach 

From the kernel on the cylinder obtained in the section Eq.( |3.36D , the partition func- 
tion on the torus can be readily obtained by sewing together the two ends of the cylin- 
der. Namely one takes the trace of the kernel on the cylinder with the appropriate measure 

Z[t, 6,] = / dUK[T, 4; U, U] = N\Vn (^) ^ T^h^f 
exp|-27rznl[fc- ^^^^~^^ ]| £ 5 + 1^ fj \{ d\J [j^ X^j 
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exp 



N 



T 



2TTn 



(3.49) 



To calculate the integral in (|3.49|) we can proceed as in the previous section, first writing the 
S{J2a ''^ct + n) as an integral, completing the square in and then Poisson resumming. One 
gets 



Z[r,9k] = N\ 



1 

2^ 



+00 



{na}=-oo 



n=0 P 



•^0 a=l \a=l I 



2-K 



(i?9exp 



Ne^Lr 
87r2 



N 



E 

\ 2 



iV 



/ N 

E 



TV 



0=1 



a) y^>a 

(3.50) 



We can now shift the integer —>■ ria — — N = Va a = 1, . . . ,N and use the sum over n n 
to complete the Gaussian integral in {}. The S{J2a^a) can be eliminated integrating in Aat. 
The integrals in Aq,, a = 1, . . . , — 1 give just a product of Kronecker deltas imposing the 
N — 1 conditions 



Ta-rN - rp(a) + rp(N) = 



a = 1, 



AT- 1 



(3.51) 



these can actually be extended to a = A^, because the is trivially satisfyed. The sum of 



the conditions (|3.51|) , for a = 1, . . . 
give rise to the determinant factor 



A^, gives r^r — rp(Ar) = 0, so that the product of Kronecker 



E(-i)^ n 



rp{a),1'a 



det 5. 



P a=l 

The SU{N) partition function on the torus then reads 

N~l +00 f 

^ exp < - 

2 



a/3 



(3.52) 



Z[r,4] = Ar!(2vr)(^-i)/2 5: ^ 

n=0 {ri,...,rjv_i}=-oo 



-27rm— 
A^ 



k- 



N 

E^° 

a=l 



det 

«/3 



'73 



exp 



e'Lr 



N 



E-«- 



a=l 



1 

N 



/ N 

E 

\a=l 



(3.53) 



The determinant of Kronecker deltas forbids any two of the integers a = 1, . . . ,N from 
taking the same value and, the expression ( p.53| ) being symmetric under the permutation of 
any r^, the integers can be ordered according to ri > . . . > = —N provided we multiply 
by a factor A^!. 

The result ( p.53|) gives the correct spectrum for 2-dimensional SU (N) Yang- Mills theories 
on the torus which, without the theta angle, has already been provided in many papers [0, |l|, ||; 
it does not agree with the one found in ref.[|l^. Zsu(n) is a sum over the representation 
of the exponential of the second Casimir, sum which is actually restricted only to those 
representations whose Young tableaux have rows with k mod[A^] boxes, where k is the discrete 
theta angle of 2-dimensional Yang-Mills theories. To compute (|3.53|) we followed the path 
integral calculation in [O generalizing it to the case of non trivial theta states. 
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4 Loop Correlators on the Cylinder 



In this section we shall compute, on the cylinder, the correlator of two Wilson loops in the 
fundamental representation. One loop is situated at the point u with / windings and and the 
other, at the point v ,with m. 



Pi 



Lm \ 



Ui,U2]U,v) 



Ui TT{U\u)}TT{U"'{v)}Pk 



(4.1) 



where Z^ coincides with the kernel on the cylinder in the k ^-sector {e.g. Eq.( |3.3(j| )). This will 
allow us to obtain also the correlator on the torus just by sewing together the two ends of the 
cylinder, as we did for the partition function in the previous section. Using the prescriptions 
introduced for the partition function on the cylinder, the path integral representation of the 
correlator (|4.1|) reads 



Pl,ra{0k,Ui,U2;U,v) 



1 



N^Zk 



-indi. 



te[o,T] 



Tr{f/'(n)}Tr{f/"(t;)}4i (f/(0), U^) 6,, f/(r), e^^/^f/a ^(t/(0))^(f/(r)) 



(4.2) 



Taking, for the time being, v > u, in order to compute ( |4.2| ) it is most convenient to perform 
first the integration in the open intervals t G (0,u), {u,v) and {v,t). This corresponds to 
the calculation of the partition function on the three cylinders with boundaries {Ui,U{u)), 
{U{u),U{y)) and {U{v),U2). After the diagonalization and the integration on the gauge 
potential in the whole interval t G [0, r], we can perform the path integral in the open intervals 
t G (0,m), (M,f) and (f , r) as we did in the previous section in the open interval t G (0,r). 
The results will then be given by Eqs.( p.2i^ ) and (|3.33| ) with the appropriate changes for the 
different lengths of the three cylinders. Then, using Eq. ( |3.35| ) for the delta-functions in ( ^.2| ) 
one can integrate on the boundary points t = 0, t = t. The result reads 



Pi, 



1 / 



(4.3) 



AT-l 



NZk \ {e^Lfu{v -u){t~v)^ 



EE 

n=0 P,P' 
N 



(^_Y^P+P' ^-2iTin[k-N{N-l)/2\/N 



J(Ai)J(A2) 

N 



+00 { ^ \ { ^ \ r ^ ^ 

E '5 E^"M E^« + ^ / X{d^''{u)d^'^{v) y: 

{na,la}=-oo \a=l / \a=l / a=l /3,7=1 

^Ev^-^l^fE 

\a=l / \a=l 



(f {v) exp 



1 ^ 



exp 



N 



2e'^L(v-u] 



E(^^ 



a=l 



N 



2e^L(T - V 



E^ 



a=l 



Af^"^ + 27r/, + 271^ 



where we have used the delta-functions in (|4.3| ) to eliminate the matrix Q appearing in ( |3.29|) . 
We can now proceed as before: write the delta-functions in the integers as in Eq.( ^.37 ), 
perform a generalized Poisson resummation in ria and la to the new integers n'^ and l'^, 



introduce the integers Tq, 



I'M ' 



■ N and Sr 



l',-Vj,-N{rN 



-N) and integrate 
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on the variables introduced to impose the constraints on the integers (say d and d') using the 
sums over n'^ and l'^ to complete the Gaussian integrals in d and d' . 



Pl,ra{0kAlA2]U,v) 

1 / 1 



(4.4) 



N-l , 

7V-1 +00 



e^Lu 



exp 



exp 



exp 



N 



1 



N 



n=0 {r„,s„}=-oo P,P' 

e^L{T - v) 



P+P'^ 



-2mn[k-N{N-l)/2-J2^ Sc]/N 



J(Ai)J(A2 



\a=l 



N 



1 



N 



Ka=l 



N 

Q = 



P(q) .P'(a) 



N / N 



N 



o=l /3,7=1 \a=l / \a=l / 

./^'^(w) + ^m^^(t;) - I E (r„^"(w) - 3„^"(t;)) - —pr-, r E (^°(^) " ^"(^))' 



Taking advantage of the delta-functions in Eq. (|4.4|) to integrate on (p^{u) and (p^{v) the 
integral in ([4 .41) becomes 



N-l N 

n difi''{u)difi"{v) E exp 

a=l /3,7=1 



N-l 



E - '^Af - ^(^a/3 - (^A^/?)) </'°(m) 



a=0 



exp 
exp 



7V-1 



L a=0 



N-l 



2e2L(t; - Uj ^^^^ 

We can now complete the square, for example in ip'^{u), using the inverse of the matrix Q 



(4.5) 



-1 

-1 N-l 



-1 \ 
-1 



V -1 



(4.6) 



-1 ... 



and then perform the Gaussian integration. At this point the integrals on the ip"'{v) just give 
products of Kronecker deltas relating the integers Sq, to the Tq, through the conditions 

Sa = ra- n{Sal3 " Sn/s) - m{6a^ - 6N-y) ■ (4.7) 

Summing on the Sa, we finally get, for the SU (N) correlator of two loops on a cylinder {v > u) 

+00 Af-l 



Pl,m{dkAlA2\U,v) 
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E ('"aAf - (r„ - IS^f^ - 'mSa-y)X2 ^°'') 
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e^LT 



exp 
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E '^a ~ 2mr^ + 2lm ( 5, 
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N 
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the case v < u can be obtained by exchanging u ^ v. For / = m = we get, as we 
should, 1. Even if in (|4.8| ) the sum over the integers Va is unrestricted, the presence of 
the determinant Z]p(~l)'^ exp(i '^a-^f forbids any couple of the to take the same 
values. Consequently, one can easily reconstruct in ( |4.8| ) the second Casimir and the sum over 
representation as in Eq. (|3.41|) . The sum will then be restricted only to those representations 
whose Young tableaux have rows with a number of boxes given by A; + / + m = 0, mod[A^] 
where k is the theta angle and /, m are the number of windings of the Wilson loops. 



4.1 Wilson and Polyakov loop correlators on the torus 



From Eq. ([4.8|) we can easily get the correlator of two Wilson (and, by duality r ^ L, Polyakov) 
loops in the fundamental representation on the torus, just by taking the trace of Eq.( [4.8D 
according to 



I n fx: A,) Pi 

•' a=l \a=l ) 



',m(6'fc, A, A;M,t;) . 



(4.9) 



Using the delta-function to eliminate the A^**^ component of A, the integrations just give a 
product of Kronecker deltas imposing the conditions 



- rp(a) + rp(^N) - tn ~ I {^a^ - Sn/b) - m {6a^f - S^^) = a = l,...,N 



(4.10) 



where the N^^ equation does not arise from the integration but can be added because is 
automatically satisfied. The sum, from 1 to A^, of these conditions gives 



rp{N) -fN- j^im + + 16ni3 + ^6^^ = , 



(4.11) 



consequently m + / must be an integer multiple of N or must be zero otherwise the delta 
imposing the conditions ( 4.10| ) will set to zero the correlator. This is nothing but the Z^v 
invariance: all the correlators that are not Z^v-invariant must vanish. The correlator (14.91) is 
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(4.12) 



Setting 171 = from Eq. ( [4.12| ) one gets the correlator for a single loop. From Zn invariance, 
this will be different from zero only if / is an integer multiple of N. 
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exp 
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(4.13) 



For a loop anti-loop correlator, m = — /, Eq.( ^4.12 ) can be rewritten in a different form to 
make clear its invariance under translations on the torus 
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(4.14) 



Eqs.( |4.1^ ) is different from the corresponding one obtained in Ref.|T^. As a matter of fact 
the restriction on the integers Tq introduced by the determinant factors, not only changes the 
spectrum of the theory but also alters the asymptotic behavior of the correlators. Consider 
for example the case of the correlator of a quark- antiquark pair, / = 1. As is known, any non- 



trivial theta-sector is unstable [0, |3^, Q in this case since the quark-antiquark pair is ener- 
getically favoured to reduce the associated background electric field. The asymptotic behavior 
of Eq. ([4.14|) can be obtained by sending r — > oo. In this limit only the lowest representation 
survives in the sums in Eq.( |4.14| ). The lowest value for the quantity I]a(ra)^ — (J^a^a)'^/^ 
can be obtained by choosing = —a (which corresponds to the singlet representation) and 
is given by the curvature of the group manifold A^(A^^ — 1)/12 as in Eq.( |3.42| ). The large r 
hmit of Eq. ([4.14| ) is 



lim Pi„i(u,i;) 
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20---N 



exp 
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(4.15) 



We then have in this sector a confining potential between the quark-antiquark pair with string 
tension e^{N'^ -^)/N. 



4.2 Results for SU{2) 

Let us now study in particular the SU{2) case on the torus. For N = 2 Eq. ( [4. 141) reads 
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Recall that the integers r + 1 for r > — 1 or — r — 3 for r < —3, in Eg. ( |4.16| ) give the number 
of boxes in the rows of the Young table in a given representation. Consequently, for k = 
only representations whose Young tableaux have an even number of boxes are present in the 
sum, for k = 1 only those with an odd number. Identifying j = (r + l)/2 for r > —1 and 
j = — (r + 3)/2 for r < —3, j = 0, 1/2, 1 ... we can rewrite the loop correlator for SU{2) as 



Pi,.i{k = 0;u,v) 
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4 J2 
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(4.17) 



and 



Pi-i{k = l]u,v) 



cosh 



e^L{v - u)l{2j + 1) 




1 l'^\v — u\ 

J J + 1 + - + ^ — ^ 



\v — u\ 



(4.18) 



Taking r — oo in the fundamental representation with / = 1, we can see that the non trivial 
theta sector has an instability which creates a repulsive potential between the external charges. 
Moreover, the prescence of fundamental charges results breaks the center symmetry that gives 
rise to the theta-sectors in the first place. For example, with a pair of fundamental charges, 
one region of the torus carries fluxes with even numbers of boxes while another carries fluxes 
with odd numbers of boxes. In general one should sum over all theta-sectors in this case to 
obtain correct results for correlators. Turning to asymptotic behaviour, using Eq. (|3.53|) it is 
easy to see that the 9 = sector produces the leading behaviour with inter-quark potential 
given by 

1 3e^ 

--logPi,_i(0;M,f) = — \v-u\ . (4.19) 
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r— >oo 
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Let us consider now the correlator for a pair of adjoint SU (2) loops. A loop in the adjoint 
representation can be taken as the modulus squared trace of the fundamental representation 
group element 

Trf/adj(n) = |Trf/(n)|'-l . (4.20) 

For SU{2), 

Trf/,dj(^i) =Trf/2(n) + l . (4.21) 



For the correlator we then have 

(Tr?7adj(^x)Trf/,dj(t;)) = 1 + {TtU\u)TtU\v: 
Using Eqs. (|4A2|) and ( CT) 
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+ exp 



e^L\v - m| (r + 2)] I - ^ exp [-e^Lr + e^L\v-u\- mn{k - 1)] I (4.23) 



n=0 



As before, only representations with an even number of boxes, r even, will contribute when 
/c = and representations with an odd number of boxes when k = 1. Taking into account the 
behavior of Zk-, Eg. ( ^.531 ), for large r, we obtain a confining behavior in the trivial theta-sector 



lim -- log (Trf/adj(u)Trt/adj(t^))o = 2e^|^^ - u\ , 
and a screening behavior in the non-trivial theta-sector. 

Jirn -y log(Tr[/adj(M)Trf/adj(f))i = j 1 - exp (-3e^L|t; - m 



(4.24) 



(4.25) 



These results are in agreement with the general discussion of ref . ||3^ • 

As a check of previous calculations, we will now evaluate these correlators for the case of 
SU{2), using the character representation. We start by considering the general case of the 
correlator of a pair of Polyakov loops in (irreducible) representations R and R' on a cylinder 

Pr,r'{Uu U2; t;) = [dVi][dV2]K[u; f/i, Vi]xr{Vi)K[v - u; Vu V2]xr'{V2)K[t - v; V2, U2] 

(4.26) 

Using the definition of the kernel K (|3.11|) and the properties of the group characters this 
expression can be reduced to a sum over representations foTV>u 

Pr,r'{Ui,U2;u,v) = (4.27) 

^ E N^}iA'XRAUi)XRsiU2) e-i<^^^[^^^(«0+(-«)c.(i..)-f(.-.)c.(«3)]} 
Ri,R2,R3 



Using this result (f4.27|) we can easily find the analogous on the torus. We can set Ui = U2 = U 
and integrate on U to obtain the pair correlator of Polyakov loops on the torus Using the 
properties of characters this reduces to 

PR,R'iu,v) = ^ Y: N^^^j,Ng',n'^w{-e'L\v-u\C2iRi)-e'LiT-\u-v\)C2iR2)) (4.28) 
^ R1R2 

From this general formula one can immediately make quantitative statements about the 
binding between pairs of loops in the case where one side of the torus becomes large (r 00). 
In the trivial theta-sector (recall this means we consider the sum of all distinct theta-sectors) 



(4.29) 



As well one can easily do the same for the correlator of a pair of Polyakov loops in the k^^ 
theta-sector 
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In the case of two Polyakov loops in the adjoint representation some calculation ||33| shows 
for — 1 > > 1, the pair correlator is 

The cases k = 1 and k = N — 1 are given by excluding the second and third terms respectively. 
The case k = recovers the topologically trivial case ( [4.29| ) in this limit. 

Let us now consider in particular the SU{2) case where the fusion numbers are known 
explicitly. We will label representations by the single non-negative integer n which is equal to 
li in terms of row variables or mi in terms of column variables in the associated Young table. 
Consequently, C2{n) = {{n + 1)^ — l)/4. The fusion numbers are 

m = ll when^=j + Z j + /-2,...,|j-/| 
■' [(j otherwise 

Bringing these facts together we have the correlator of pair of Polyakov loops, one in 
representation n and the other in m, separated by distance a. \v — u\ 

PnAu^v) = (4.33) 



oo 2ri I 2, mill i,m min jr,s 

z 



j,l=0 r=0 s=0 

The first delta function serves to enforce the condition that m + n must be even. This is a 
particular example of the general fact in SU (N) that in order for the correlator of any system 
of loops to be non-vanishing the total charge of the loops must be vanishing mod A^. This 
restriction ensures that the system contains a charge singlet. 

Two special cases that are of interest are the pair correlators of fundamental and adjoint 
Polyakov loops. For the fundamental case n = m = 1 and 

v) = lf: e-^[('+i)^-i][ e-^(2'+3) + e- ^'''^4'"-"'' (^^+^)] (4.34) 

Zj 



1=0 



It can be checked that this result coincides with the sum of Eqs. ( [1.17| ) ( [4.18|) when the 



different zero point energy and normalization of Z in the two approaches, are taken into 
account. Likewise for the correlator of a pair of adjoint loops m = n = 2 

P2,2{U,V) = ^ ^-2e^Lir-\v-u\) ^435) 

Zj ^ 

+ ^ Q-^-T^imf~i]^ Q-eU\v-u\{l+2) ^ ^ ^ g-e2L(r-|^-«|)(Z+2)| j 
1=1 J 

Now we turn the issue of non-trivial ^-sectors in the correlators. As usual this is carried out 
by including a projection operator. For SU{2) this prescription amounts to restricting the sum 
over representations to those with either odd or even numbers of boxes in the corresponding 
Young table. It should be noted that this can be done in a consistent manner only if each of 
the charges in the system has vanishing 2-ality (ie. has 2n boxes). 
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For the correlator of a pair of adjoint Polyakov loops the two cases are k = where one 
would sum over even representations in ( [4 .331) 

00 2 \ 

i=2,4,--- / 

and the case k = 1 where the sum is over odd representations 

P2,2iev,u,v) = ^ E e-^[('+i)'~i][e-^'^l«-^l('+2) + i+ e-^'^(^-l^-"l)('+2)] (4.38) 

^1 /=1,3,- 

These results coincide with the results of the functional integral technique, Eg. ( [4.23| ). 

5 Discussion 

We have demonstrated the utility of the diagonalization method for doing practical com- 
putations in 2D Yang-Mills theory. This method is alternative to the well-known charac- 
ter expansion technique. The diagonalization method is useful for computing correlators of 
higher winding loop operators. Knowing the formulae for higher representations in terms 
of higher winding loops (via characters), and the formulae for correlators in the character 
representation, it should be possible to obtain information about the fusion rules for higher 
representations of SU(N) groups (and in fact, by suitable generalization of our work, arbitrary 
Lie groups) from our formulae. We have not pursued this problem in the present work. 
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